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Abstract 



23 ] in two di- 



Building on the recent solution for the spectrum of the non-commutative well 
mensions, the thermodynamics that follows from it is computed. In particular the focus is put on 
an ideal fermion gas confined to such a well. At low densities the thermodynamics is the same as 
for the commutative gas. However, at high densities the thermodynamics deviate strongly from 
the commutative gas due to the implied excluded area resulting from the non-commutativity. In 
particular there are extremal macroscopic states, characterized by area, number of particles and 
angular momentum, that correspond to a single microscopic state and thus have vanishing en- 
tropy. When the system size and excluded area are comparable, thermodynamic quantities, such 
as entropy, exhibit non-extensive features. 

PACS numbers: ll.10.Nx 
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There seems to be growing consensus that our notion of space-time has to be drasti- 
cally revised in a consistent formulation of quantum mechanics and gravity 

HQ. 

One 

possible generalization, suggested by string theory that has attracted much interest re- 
cently is that of non-commutative space-time Despite a number of investigations into 
the possible physical consequences of non-commutativity in quantum mechanics and quan- 



tum mechanical many-body systems [5|, [6, |7|, |8j, quantum electrodynamics [9|,llQ|,lll(, the 
standard model [ijl and cosmology [l3,[l4]), our understanding of the physical implications 
of non-commutativity is still in its infancy. In particular our understanding of the physics 
is cluttered by the implementation of the Poincare symmetry. One either has to accept 
explicit breaking of the Poincare symmetry, with its possible physical ramifications in e.g. 



high energy cosmic radiation 
the Poincare group 
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or one can restore the Poincare invariance by twisting 
In the latter case it is also necessary to twist the 



statistics, but there is a controversy on whether this should be done in the sense of 
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2Q|, 



and it 



which leads to modified statistical correlations [20|, or the braided sense of 
seems difficult to distinguish between these two possibilities. Indeed, the physical effects 
of non-commutative space-time and modified statistical correlations arising from twisted 
statistics seem to be irrevocably intertwined. Even in the case of a non-interacting, but 
confined non-commutative gas, the physical consequences of non-commutativity have not 
been investigated in detail due to the problem of giving precise meaning to boundaries in a 
non- commutative space. Against this backdrop it is highly desirable to find systems where 
the physical effects of non-commutativity can be investigated exactly and unambigously. 
Recently precise meaning was given to a particle confined to a two dimensional non- 



commutative infinite well 



231 ] . the spectrum was computed exactly and it was shown that 



the results reduce to the conventional results in the commutative limit. The result, summa- 
rized below, is that at low energies there is very little distinction between the commutative 
and non-commutative world, if the non-commutative parameter is small compared to any 
other length scales squared in the system, which will be the case if the non- commutative 
parameter is of the order of the Planck length squared. Thus one expects that at normal den- 
sities (interparticle spacings) and temperatures (thermal length) non-commutativity should 
be virtually undetectable. At high densities or temperatures this situation may, however, 
change. This motivated us to study the thermodynamics of such systems at extreme densi- 
ties and pressures, which is what we report on here. The availability of exact results for the 
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spectrum enabled us to study these limits reliably, albeit numerically. Furthermore in two 
dimensions the commutator [x,y] is already a 5*0(2) invariant ([x',y f ] = [x,y]) as one can 
easily check and thus the non-commutative relation [x, y] = i9, where 9 is considered to be a 
5*0(2) invariant, does not break the rotational symmetry. In two dimensions it is therefore 
not necessary to restore the rotational symmetry by twisting the rotation group. This is in 
contrast to the higher dimensional case where the non-commutative relation = in- 

breaks the rotational (Poincare) invariance if the 6^ are considered to be fundamental invari- 
ants which are the same for all observers and a twist of the rotational symmetry is required 
to restore the rotational invariance. The complication of a twisted rotational symmetry, 
twisted statistics and possible modified statistical correlations arising from it does therefore 
not occur in two dimensions, which opens up the possibility of studying the physical effects 
of non-commutativity in an unambiguous setting. There is a possible exception to this in the 



case of a semi-direct product of a twisted gauge symmetry and rotational symmetry [24j |. but 
this does not apply to our current investigation. Furthermore it should be pointed out that, 
apart from the issue of the restoration of rotational symmetry and the associated question 
of statistics, the current analysis can be extended straightforwardly to higher dimensions. 
Indeed, in three dimensions one can, by an appropriate change of coordinates, always re- 
duce the non-commutativity to two non- commutative and one commutative coordinate. As 
expected this breaks rotational invariance, but if one does not insist on its restoration, the 
analysis given below extends with minor modifications to the higher dimensional case. This, 
and the associated issue of the restoration of rotational symmetry, will, however, be the 
topic of future investigations. 



We start by summarizing the key results of 23( which we use here. The spectrum of a 
particle with mass mo confined to a disc of radius R 2 = 6(2M + 1), with [x, y] — id (9 > 0), 
in a non-commutative space is given by 



Ok 

LZ +1 (-^1=0, m > 



L l M ] -\m-< ('—) = °> -M<m<0, (1) 

h 2 k 2 
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Here M is an arbitrary non-negative integer that fixes the radius of the disc, m is the 
angular momentum and an associated Laguerre polynomial. Eigenf unctions for m < —M 




FIG. 1: Spectrum of a particle in a non-commutative well for M = 10 in units of Eq = as a 
function of angular momentum. 

vanish identically and the spectrum truncates at angular momentum m = —M. A crucial 
observation, underlying many of the features discussed below, is that each positive angular 
momentum sector has exactly M + 1 excitations and each negative angular momentum 
sector, — M < m < 0, exactly M — \m\ + 1 excitations. For later reference we give a 
graphical summary of these results in Fig. 1. 

The spectrum for the corresponding commutative system is given by 

J m (kR) =0, m = 0,±l,±2... 

E = (2) 
2m V ; 

with J m a Bessel function. 

To compute the spectrum and thermodynamics of the system, we must therefore be 
able to compute the zeros of the Laguerre polynomials and Bessel functions. While efficient 
algorithms to determine the zeros of the Bessel functions are available, the most efficient and 
accurate determination of the zeros of the Laguerre polynomials is to find the eigenvalues of 
the following N x N hermitian matrix, constructed such that its eigenvalues yield exactly 
the zeros x r ^ m (1 < r < N) of the Laguerre polynomial U£{x): 



H itj (N, m) = (2i - 1 + m) 5 itj + y/i(i + m)S j>i+1 + y/j(j + m)5 j ^ 1 , i, j = 1, 2, . . . N. (3) 

Indeed, the characteristic polynomial in x of this matrix obeys a three term recurrence 
relation which is that of the Laguerre polynomial L™(x) up to a normalisation factor, 
det (H(N, m) - xl) = N\ L%(x). 



From Fig. 1 it should be clear that for certain choices of total number of particles and total 
angular momentum, the single particle states must be filled in a unique way (for modelling 
purposes in the present paper, these particles are taken to obey a Fermi statistics with a 
single spin state). Indeed, let us consider the situation where all the single particle states 
from angular momentum — M up to angular momentum m are completely filled. The total 
number of particles is 

m+M 



N„ 



^2 {k + 1) = (1 + m + M) ( 1 + - (m + M) ) , —M < m < 0, 
k=o ^ ' 

u-\ ^ 

^2(k + 1) + (m + 1)(M+ 1) = - (1 + M) (2 + 2m + M) , m > 0. (4) 



fc=0 

The total angular momentum of the state is 

m+M 
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L m = V* (k + 1)0 - M) = -{2m-M) (1 + m + M) (2 + m + M) , -M < m < 0, 

fc=0 

M— 1 m 

= J^(A; + 1)(A;-M) + ^(M + 1)A; 

fc=0 fe=0 

= -(1 + M) (3 m (l+m)-M (2 + M)), m > 0. (5) 

It is clear that any other distribution of the number of particles (T4]) must lead to a higher 
angular momentum than ([5]), which is thus the minimum total angular momentum for this 
number of particles and system size M. This implies that if we fix, for system size M, the 
total angular momentum on ([5]) and the number of particles on (T4j) for some m, there is only 
one possible distribution of particles. The entropy for this macroscopic state thus vanishes. 
Furthermore it is clear that increasing the particle number, while keeping the system size 
fixed, or decreasing the system size while keeping the particle number fixed, must yield 
angular momentum values larger than L m . Thus, if the total angular momentum, which is 
a conserved quantum number, is fixed on L m , we expect these states to be incompressible 
and the pressure to diverge as either the particle number or system size approach the critical 
values from below or above, respectively. Alternatively we expect the pressure to diverge 



as the particle density approaches the critical value p m = nd (2M+i) f rom below. The total 
energy of these states can also be evaluated quite easily and in units of E = is given 
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by, 



„ m+M 

Hi, 



=^ = V trH(k + l,M- k) 



k=0 



= ( ! + \ f) ( 1 - + \/) ( 1 + ^ (m + M) ) . -1/ < //; < 0. 



M-l 



= ^tr#(£; + l,M-£;) + ^tr#(M + l,£;) 

fc=0 fe=0 

= i (1 + M) (M (1 + M) + (1 + m) (2 + m + 2 M) ) , m > 0. (6) 

For the extremal cases above we have considered the micro canonical ensemble, i.e., 
fixed area, particle number, total angular momentum and energy. It is, however, only in 
these extreme cases that that the computation of the micro canonical partition function and 
thermodynamics is trivial. For general particle numbers and angular momenta this becomes 
very difficult and to proceed in these cases it is more appropriate to switch to the grand 
canonical ensemble in which the average particle number, average energy and average total 
angular momentum are fixed. The reason for fixing the average total angular momentum 
is to ensure that when the commutative and non-commutative systems are compared, it is 
done for the same physical situation. It should be clear from the spectra of the commutative 
and non-commutative systems that fixing only the particle number and temperature may 
lead to vastly different average total angular momenta. 

The appropriate density matrix for this ensemble that acts on Fermi Fock space is 

£ = tre -/^-^-^). ( 7 ) 

Here the Hamiltonian, H, particle number operator, N, and total angular momentum, L, 
are given in second quantized form by 

r,m 

N = ^X m a r , m , (8) 

r,m 

L = J2 hma t,m a r,m, 

r,m 

where aj m , a r>m creates (annihilates) Fermi particles with energy E r>m and angular momen- 
tum m. Furthermore (3 = 1/kT, p is the chemical potential and uj a Lagrange multiplier 



that enforces the constraint of fixed total angular momentum, which, for rigid bodies, has 
the physical meaning of the angular velocity. The partition function, £, which involves a 
trace over full Fermi Fock space, is explicitly given by 

£ = ^ Te -l3{H-nN-uiL) = TT (l + Q-^^m-^-nwm)^ _ ^ 

r,m 

The central thermodynamic quantity from which all thermodynamics derive, is the q- 
potential which is the logarithm of the partition function 

q{A, T, fi, w) = log £ = lo S i 1 + e-^r^-v-^m)^ _ ^ 

r,m 

The q-potential depends on the temperature, chemical potential, angular velocity and area, 
A, of the system. The dependence on area enters through the single particle energies that 
depend on it. 



To find the thermodynamic interpretation of the q-potential one can proceed as in 26] 
p. 101 with slight modification due to the constraint on angular momentum. For completeness 
we briefly outline the argument. Taking the differential of the q-potential one has 

J TT 

dq = -d(3 (E — jjN — uoL) - (3(—}dA + fiNd/i + (3Lduo (11) 
where E, N, L and (4j) denote the average energy, particle number, angular momentum 



and average of respectively. This can be rewritten as 



d H 

d{q + f3 {E - fiN - ujL)) =(3\dE- fidN - udL - (~^)dA ) . (12) 



Comparing the right hand side with the first law for a rotating body 25| p. 74 



TdS = dE + PdA — fidN — udL, (13) 
one concludes that P = — and 

1 Q 

d(q + (3{E-fiN-ujL)) = —, (14) 

K 

or, 

9 

q = - - /3(E- fiN -cuL). (15) 

I — I 

Keeping in mind that the Gibbs free energy for a rotating body is given by [25] p. 74 

G = E — uL — TS + PA = /J.N, (16) 



we have 

"=w- < 17 > 

The rest of the thermodynamics follow easily with the quantities most important to our 
analysis given by 

N = kTp-, 

L = (18) 
dT 

We are now in a position to systematically compute the thermodynamics of the commu- 
tative and non-commutative Fermi gas. To do this, and for comparison, it is convenient 
to first rewrite the expression for the q-potential in terms of dimensionless quantities. We 
start with the non-commutative case. Let x r ^ m denote the zeros of the Laguerre polynomi- 
als Lm +1 (x r>m ) = 0. From flTJ we then have E r>m = - ''q^™ ■ Introducing the energy scale 
Eq = and the dimensionless parameters (3 = E (3, jl = -^h and u = |r, the q-potential 
for the non-commutative well can explicitly be written as 

oo 

q(Mj,ji,u>)= ^log(l + e-^--^). (19) 

m=—M r 

The dependence on the area of the system is reflected by the dependence on M, keeping in 
mind the relation A = nR 2 = tt9(2M + 1). 

The commutative case is treated in a similar fashion. We denote by j rjm the zeros of the 
Bessel functions J m (j r ,m) = 0- From (j2J) we then have E r>m = 2m r Q ^ 2 ■ Using the relation 
R 2 = 8(2M + 1), the q-potential for the commutative well can be expressed as 

q(M, 0, fi,Cu) = J2 lo S f 1 + e~^^2-^ m) \ . (20) 

r,m ^ ' 

In terms of these dimensionless variables the thermodynamic relations (fT8|) read 

n = 

L 1 dq 

S ~ dq . „. 

= q-P-1. 21 
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FIG. 2: Average number of particles as a function of the chemical potential at temperature (3 = 1, 
system size M = 10 and average angular momentum L = 0. The solid line is the non-commutative 
case and the dashed line the commutative case. 

We start our analysis by noting a fundamental difference between the commutative and 
non-commutative systems. In the discussion of the extremal states, we have observed that 
there is a minimum total angular momentum for a given system size and particle number. 
Thus, for a fixed angular momentum and system size, M, the number of particles in the 
non-commutative system can not be made arbitrarily large. When all the negative angular 
momentum states are filled, the angular momentum must necessarily increase with particle 
number until the minimum total angular momentum will exceed the given total angular mo- 
mentum value. This is shown in Fig. 2 where the number of particles is shown as a function 
of the chemical potential for fixed angular momentum L = 0. One observes that the particle 
number saturates at a maximum value. This implies that the non-commutative system has, 
for a given angular momentum and system size, a maximal density, while no such limit exists 
in the commutative system. This cut-off in the density for the non- commutative system is 
a reflection of the implied excluded area resulting from non-commutativity. The maximal 
density clearly has a temperature dependence, which turns out to be very weak at the tem- 
peratures we consider and we shall neglect it in what follows. Note that, not surprisingly, 
the results for the commutative and non-commutative cases agree at low densities, but that 
they start diverging strongly at high densities. 

The maximal particle density as a function of system size can easily be computed and is 
depicted in Fig. 3. Here we show the dimensionless particle density p defined by p = = 
tt0(2M+i) = tR>- sna ded region is forbidden, i.e., each system size has a maximal density 
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FIG. 3: Maximal dimensionless density as function of system size at j3 = 1 and L = 0. The 
unshaded region is physical. 

associated with it. Conversely, a fixed density can only be accommodated by systems with 
a size greater than a minimal size. We can also find an analytic estimate for the curve 
shown in Fig. 3. From ([5]) we can easily establish the value of the angular momentum, 
m, up to which all levels must be filled in order to yield L = 0. In general this is not an 
extremal situation, but it is simple to see that this only introduces correction of higher order 
in 1/M. Thus, to leading order in M we find m « Substituting this in (j3J) we find for 
the dimensionless density p « | ^1 + -^j M = 0.539M. The corresponding slope in Fig. 3 
is 0.536. Although our analytic estimate of this slope is done at zero temperature, while 
the computation in Fig. 3 is done at $ — 1, we still expect good agreement as (3 = 1 is 
of the same order or less than the dimensionless single particle level spacing, i.e., we are 
approximately at zero temperature. This, combined with the weak temperature dependence 
of the maximal density, lead us to expect that this temperature difference will have a small 
effect, as is indeed the case. 

Next we consider what happens in the non- commutative system as the maximal density 
is approached. As remarked before, we expect on purely physical grounds an incompressible 
behaviour with diverging pressure. In Fig. 4 we compute the pressure as a function of particle 
density for a fixed system size, i.e, the pressure as a function of the number of particles. 
We observe that for the non- commutative system the pressure does indeed diverge at the 
maximal density, also shown on the graph, while the pressure for the commutative system 
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FIG. 4: Dimensionless pressure as a function of density for the non-commutative (solid line) and 
commutative (dashed line) systems at two temperatures (3 = 0.2 and (3 = 1, system size M = 10 
and angular momentum L = 0. The critical density at (3 = 1, where the non-commutative pressure 
diverges, is indicated by the vertical dashed line. As remarked in the text, the critical density has 
a weak temperature dependence and its values virtually coincide for the two temperatures shown, 
thus we neglect this difference here. 

shows no anomalous behaviour. The pressure shown is a dimensionless pressure denned by 

n9P q 



P 



(22) 



E p(2M+l) 

We note that at low densities the commutative and non- commutative equations of state 
coincide, while they deviate strongly at high densities. 

Let us consider the dependence of the pressure on the system size. In Fig. 5 the dimen- 
sionless pressure for the commutative and non-commutative systems are shown as a function 
of system size for N = 150 particles, temperature (3 = 1 and angular momentum L = 0. 
At large system sizes the pressures agree, but they deviate strongly at smaller system size. 
At system size M = 10 the non- commutative pressure diverges, while this only happens at 
system size zero for the commutative case. As we are at low temperatures here (/3 is of the 
same order or less than the dimensionless single particle level spacing), we are in the degen- 
erate limit. For L = it is easy to see that the behaviour of the commutative gas should be 
the same as an ordinary free Fermi gas and one can quite easily estimate the behaviour of 
the dimensionless pressure (in two dimensions) in the degenerate limit to be P oc (2M+ 1)~ 2 
(see 26] p. 215). This is shown by the solid line in Fig. 5. Naively one might expect the 



behaviour of the non-commutative gas to be simply that of a gas in which particles occupy 
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FIG. 5: Dimensionless pressure as a function of system size M for the non-commutative (stars) 
and commutative (box) systems at temperature $ = 1, particle number N = 150 and angular 
momentum L = 0. The minimal system size at which the pressure diverges is indicated by the 
vertical dashed line. The solid line shows the degenerate pressure for the commutative system 
given by P = 23 685/(2M + l) 2 . The proportionality constant was determined by fitting the data 
point at M = 0. 

a finite size that gives rise to an excluded volume, e.g., the van der Waals gas. However, it 
turns out not to be the case. The behaviour of the pressure in the non-commutative case 
is not simply P oc (2M — 2M )~ 2 with M denned by the minimal area. As is clear from 
Fig. 5, the pressure still behaves as P oc (2M + 1)~ 2 at large system sizes, but it deviates 
from this behaviour at system sizes close to the minimal system size. One can, of course, 
fit curves through the points in Fig. 5, but this does not really lead to new insight and we 
refrain from doing this. It would be much more helpful if a simple analytic understanding, 
as in the commutative case, of this behaviour can be developed, but we have not been able 
to do it thus far. 

The entropy of the commutative and non-commutative systems also reflects the funda- 
mental differences at high densities. Fig. 6 shows the entropy as a function of density for a 
system size M — 10 and angular momentum L = 0. One observes that the entropies of the 
two systems coincide at low densities, but that they diverge strongly at high densities and in 
particular that the entropy of the non-commutative system decreases. This behaviour can 
easily be understood by looking at the non-commutative spectrum from which it is clear that 
the density of states available to the system at high densities and fixed angular momentum 
must decrease. Indeed, for the extremal cases discussed earlier the entropy vanishes at the 
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FIG. 6: Entropy in units of k as a function of density for the non-commutative (solid line) and 
commutative (dashed line) systems at temperature (3=1, system size M = 10 and angular 
momentum L = 0. The critical density at (3 = 1 is indicated by the vertical dashed line. 

critical density. 

Fig. 7 shows the entropy as a function of system size for particle number N = 150, 
temperature (3 = 1 and L = 0. As in the case of the pressure, the behaviour of the entropy 
of the commutative system can again be understood easily as the degenerate limit of a free 
two dimensional Fermi gas. One can then easily estimate, for a fixed number of particles, 
(see j3] p. 215) that f oc (2M + 1), which is clearly the behaviour exhibited in Fig. 7 as 
indicated by the dashed line. The deviations at very small system size are due to finite size 
corrections. Note that the linear dependency of the entropy on system size is a particular 
feature of two dimensions and not generic. In contrast, the behaviour of the entropy of 
the non-commutative system is vastly different and tends to zero at the minimal system 
size, indicated by the vertical dashed line. Note that the minimal system size can be of 
macroscopic scale if the density is high. At large system sizes (low density as the particle 
number is fixed) the dependency on system size tends to that of the commutative system. 
Again we lack a simple analytic understanding of this behaviour. 

At fixed density and large system sizes, one expects intensive quantities, such as pressure 
and chemical potential, to be independent of system size, while extensive quantities such 
as entropy and energy should scale linearly with system size. We now proceed to show 
that this expectation fails in the non-commutative system, even at large system sizes if the 
density is high enough. The deviation from extensive behaviour is again due to the presence 
of a minimal system size, arising from the excluded area implied by non-commutativity. 
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FIG. 7: Entropy in units of A; as a function of system size M for the non-commutative (stars) 
and commutative (box) systems at temperature $ = 1, particle number N = 150 and angular 
momentum L = 0. The dashed line depicts the behaviour of the entropy for the commutative 
system as derived from the degenerate limit of the free Fermi gas. The vertical dashed line shows 
the minimal system size. 

Fig. 8a shows the pressure for the non-commutative system as a function of system size 
M at two densities p — 5.1 and p = 4.037, temperature (5 — 1 and angular momentum 
L = 0. As one would expect the pressure increases sharply at the minimal system size, 
indicated by the dashed vertical lines, for these densities. Here we do not compare with 
the commutative system as we are at fairly large densities where the commutative and non- 
commutative results may differ significantly and nothing can really be learned by such a 
comparison. However, we do want to point out that the pressure in the commutative system 
only exhibits a weak dependency on the size of the system as shown in Fig. 8b for the same 
parameter values as in Fig. 8a (note the different scales on the vertical axis in Figs. 8a 
and 8b). As mentioned before, this can be expected as pressure is an intensive quantity 
and should therefore be independent of system size when temperature and density are held 
fixed. Indeed, in the degenerate limit the behaviour of the pressure as a function of density 
can easily be established to be P oc p 2 . This is of course only true at large system size 
when finite size corrections are unimportant and there may be corrections at small system 
size as one indeed observes (see Fig. 8b). In contrast, in the non-commutative case the 
sharp increase in pressure always occurs at the minimal system size, which may be large 
(see Fig. 2), if the density is high. Thus this deviation is not a finite size correction, but 
due to the excluded area resulting from the non- commutative nature of the system. Far 
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enough above the minimal system size, where the excluded area is unimportant, the non- 
commutative system starts to behave like the commutative one and the pressure exhibits 
the normal intensive scaling property. 
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FIG. 8: Pressure for the non-commutative (a) and commutative (b) systems at densities p = 4.037 
and p = 5.1, temperature (3=1 and angular momentum L = 0. Note the difference in the scales 
of the vertical axis in (a) and (b). 



The entropy of the non-commutative system exhibits similar scaling properties, shown in 
Fig. 9a, for densities p — 5.1 and p = 4.037, temperature (3 = 1 and angular momentum 
L = 0. There is clearly a strong deviation from extensive scaling at the minimal system size, 
while extensive scaling is recovered at sizes well above the minimal system size. In contrast, 
the entropy of the commutative system scales extensively for virtually all system sizes as 
indicated in Fig. 9b for a density p — 5.1 and all other parameters the same as in Fig. 9a. 
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FIG. 9: Entropy for the non-commutative (a) and commutative (b) systems at densities p = 4.037 
and p = 5.1, temperature (3=1 and angular momentum L = 0. The dashed lines are only intended 
to guide the eye and to emphasize the deviation from extensive behaviour for the non-commutative 
system. 
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